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The production of charm and beauty quarks in ep interactions has been mea- 
^ ! sured with the ZEUS detector at HERA for squared four-momentum exchange 

^ ! > 20GeV^, using an integrated luminosity of 126 pb~^. Charm and beauty 

quarks were identified through their decays into muons. Differential cross sections 
were measured for muon transverse momenta pi^ > 1.5 GeV and pseudorapidities 
— 1.6 < 7]^ < 2.3, as a function of p^, rj^, and Bjorken x. The charm and 
beauty contributions to the proton structure function F2 were also extracted. 
The results agree with previous measurements based on independent techniques 
and are well described by QCD predictions. 
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1 Introduction 



The measurement of charm and beauty production in deep inelastic scattering (DIS) 
provides a stringent test of quantum chromodynamics (QCD) since the large quark masses 
provide hard scales that make perturbative calculations applicable. At leading order, 
heavy quarks (HQs) are produced in DIS via boson-gluon fusion (BGF) {'-j*g — ^ qq). A 
precise measurement of HQ production in DIS therefore provides a direct constraint on 
the gluon parton density function (PDF) of the proton. 

Charm production in DIS at HERA has been measured previously using reconstructed 
charmed mesons [H [2] or inclusively by exploiting the long lifetime of charmed hadrons [5] . 
Beauty production in DIS has been studied in events with muons and jets [H [5] and from 
lifetime information [5]. The existing data are generally in good agreement with next- 
to-leading-order (NLO) QCD predictions. The largest differences were observed in the 
muon analyses [H [S] where the measured beauty cross section was about two standard 
deviations above the theoretical expectation. 

In this paper, a simultaneous measurement of beauty and charm production using semi- 
leptonic (SL) decays into muons is presented. The fractions of muons originating from 
charm, beauty and light flavours (LF) were extracted by exploiting three discriminating 
variables: the muon impact parameter, the muon momentum component transverse to 
the associated jet axis and the missing transverse momentum, which is sensitive to the 
neutrino from SL decays. 

The analysis focused on data with large squared four-momentum exchange at the electron 
vertex, Q^, where charm measurements based on muons are competitive with those based 
on identified charmed mesons. 

The cross sections for muons from charm and beauty decays were measured for > 
20 GeV^, muon transverse momenta pi^ > 1.5 GeV and pseudorapiditie^ —1.6 < rj'^ < 2.3 
as a function of p^, r/^, Q^, and of the Bjorken scaling variable x [B] and compared to 
QCD predictions. The muon cross sections, measured in bins of x and Q^, were used 
to extract the heavy quark contributions to the proton structure function F2 which were 
compared to previous results and to QCD predictions. 

The data used in this analysis were collected with the ZEUS detector in the 2005 running 
period during which HERA collided electrons with energy Ee = 27.5 GeV with protons 
with Ep = 920 GeV corresponding to a centre-of-mass energy a/s = 318 GeV. The corre- 
sponding integrated luminosity was C = 126.0 ± 3.3 pb~^. 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing towards the 
centre of HERA. The pseudorapidity is defined as ?/ = — hi (tan |). where 6 is the polar angle. 
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2 Theoretical predictions 



Heavy quark production in DIS has been calculated at next-to- leading order (0(q!^)) 
in the so-called fixed fiavour number scheme (FFNS) in which only light fiavours are 
present in the proton and heavy quarks are produced in the interaction [7j. The results 
of this analysis have been compared to NLO calculations performed with the HvQDis 
program [Hj [9]. The renormalisation and factorisation scales were set to /i^ = fi^p = 
+ 4m^ and the quark masses to rUc = 1.5 GeV and = 4.75 GeV. The PDFs were 
obtained by repeating the ZEUS-S [10] PDF fit in the FFNS with quark masses set to 
the same values as in the HvQDis calculation. 

To calculate muon observables, the partonic results were interfaced to a model of HQ 
fragmentation into weakly decaying heavy hadrons and of the decay of heavy hadrons into 
muons. The hadron momentum was obtained by scaling the quark momentum according 
to the fragmentation function of Peterson et al. [11] with the parameter ec = 0.055 for 
charm and et = 0.0035 for beauty. This choice of Cc corresponds to Cc = 0.035 for D* 
mesons [12] since kinematic considerations [13] and direct measurements [T^ show that, 
on average, the momentum of the weakly decaying hadrons is ~ 5% lower than that of 
D* mesons. 

The semileptonic decay spectrum for charm was taken from a recent CLEO measure- 
ment [15]. The decay spectrum for beauty hadrons was taken from the Pythia [16] 
Monte Carlo (MC), mixing direct SL decays and cascade decays through charm according 
to the measured branching ratios [17]. It was checked that the MC described BELLE and 
BABAR data [18] well. The branching ratios were set to B{c ^ fi) = 0.096 ± 0.004 and 
B{b ^ fi) = 0.209 ± 0.004 [17]. 

The uncertainty on the theoretical predictions was evaluated by independently varying 
fiR and fip by a factor two; by varying the HQ masses simultaneously to {rricmi)) = 
(1.3, 4.5), (1.7, 5.0) GeV in the calculation and in the PDF fit; by varying the proton 
PDFs by their experimental uncertainty and by varying the fragmentation parameters 
within 0.04 < ec < 0.12 (corresponding to 0.025 < ec < 0.085 for D* mesons [H]) and 
0.0015 < e;, < 0.0055. As a further check, the fragmentation was performed by scaling the 
sum of the energy and the momentum parallel to the HQ direction, E+p\\, rather than the 
HQ momentum. The total theoretical uncertainty was obtained by adding in quadrature 
the effects of each variation. In the beauty case, the total uncertainty is dominated by 
the variation of fiji and of the mass while for charm the variation of ec also gives a large 
contribution. 

The calculations of F2'^ and in the FFNS were performed using HvQDiS and cross 
checked with the QCD evolution code [20] used in the ZEUS PDF fit. 
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3 Monte Carlo samples 



Charm and beauty MC samples were generated using Rapgap 3.00 ^21j to simulate 
the leading order BGF process. Parton shower techniques were used to simulate higher 
order QCD effects. Higher order QED effects were included through Heracles 4.6 fF2\ . 
The CTEQ5L [23] PDFs were used and the HQ masses were set to rric = 1.5 GeV and 
mfe = 4.75 GeV. 

Light ffavour MC events were extracted from an inclusive DIS sample generated with 
DjangoH 1.3 [21] which is interfaced to Lepto 6.5 [25] to simulate the hadronic final 
state with the matrix element plus parton shower (meps) model and to HERACLES 4.6 
to include electroweak radiative corrections. The CTEQ5D [23] parton density was used. 

Inelastic J/ip production was simulated with Cascade [26] since that model generally 
describes the DIS data of a previous publication [27] . 

The above samples corresponded to at least five times the luminosity of the data. A 
smaller light quark sample was generated with Rapgap and mixed with the heavy quark 
Rapgap samples for the study of the inclusive DIS control sample (Section [6]). 

Fragmentation and particle decays were simulated using the Jetset/Pythia model [281 
[TB] . The lepton energy spectrum from charm decays was reweighted to agree with CLEO 
data [15]. The MC events were passed through a full simulation of the ZEUS detector 
based on Geant 3.21 [29]. They were then subjected to the same trigger criteria and 
reconstructed with the same programs as used for the data. 

4 Experimental set-up 

A detailed description of the ZEUS detector can be found elsewhere [30]. A brief outline 
of the components that were most relevant for this analysis is given below. 

Charged particles were tracked in the silicon microvertex detector (MVD) [31] and in the 
central tracking detector (CTD) [32], which operated in a magnetic field of 1.43 T provided 
by a thin superconducting solenoid. The MVD consisted of a barrel (BMVD) and a 
forward (FMVD) section with three cylindrical layers and four vertical planes of single- 
sided silicon detectors, respectively. The CTD consisted of 72 cylindrical drift chamber 
layers, organised in 9 superlayers covering the polar- angle region 15° < 6* < 164°. After 
alignment, the single-hit resolution of the BMVD was 25 /im and the impact parameter 
resolution of the CTD-BMVD system for high- momentum tracks was 100 fim. 

The high-resolution uranium-scintillator calorimeter (CAL) [33] consisted of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part was 



3 



subdivided transversely into towers and longitudinally into one electromagnetic section 
and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections. Under test- 
beam conditions, the CAL single-particle relative energy resolutions were a{E)/E = 
0.18/ Ve for leptons and a{E)/E = 0.35/ Ve for hadrons, with E in GeV. The energy of 
electrons hitting the RCAL was corrected for the presence of dead material using the rear 
presampler detector (PRES) [M] and the small angle rear tracking detector (SRTD) |35j . 

The muon system consisted of rear/barrel (R/BMUON) [3S] and forward (FMUON) [5U] 
tracking detectors. The B/RMUON consisted of limited-streamer (LS) tube chambers 
placed behind the BCAL (RCAL), inside and outside a magnetised iron yoke surrounding 
the CAL. The barrel and rear muon chambers cover polar angles from 34° to 135° and 
from 135° to 171°, respectively. The FMUON consisted of six trigger planes of LS tubes 
and four planes of drift chambers covering the angular region from 5° to 32°. The muon 
system exploited the magnetic field of the iron yoke and, in the forward direction, of two 
iron toroids magnetised to ~ 1.6 T to provide an independent measurement of the muon 
momentum. 

The luminosity was measured using the Bethe-Heitler reaction ep e^p with the lumi- 
nosity detector which consisted of two independent systems, a photon calorimeter [37] 
and a magnetic spectrometer 



5 Event reconstruction and selection 

A three-level trigger was used to select events online [301 ES]- DIS events were selected 
by requiring a scattered electron in the CAL. 

A scattered electron with energy E'^ > 8 GeV was required offline. The primary vertex 
had to be within ±30 cm in Z from the nominal interaction point. 

Muons were reconstructed by matching a CTD+MVD track to a track segment in the 
inner or outer B /RMUON chambers or to an FMUON track crossing at least four FMUON 
planes. This B/RMUON selection was looser than in some previous analyses, which 
required the muons to reach the external chambers [H HQ] , allowing a lower threshold for 
the muon transverse momentum. 

The central track associated to a B/RMUON candidate was required to pass at least 
three CTD superlayers and to have at least four hits in the MVD to allow a good impact 
parameter measurement. The tracks associated to FMUON candidates were required to 
pass at least one CTD superlayer, corresponding to at least four degrees of freedom in 
the track fit. 
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Muons were accepted in the kinematic region defined by 

pf^ > 1.5GeV,-1.6 < r]'' < 2.3. 

The hadronic system (including the muon) was reconstructed from energy flow objects 
(EFOs) jH] that combine the information from calorimetry and tracking, corrected for 
energy loss in the dead material. The EFOs were corrected using the measured momenta 
of identified muons [i0lH2]. A reconstructed four-momentum {p^x ■> Py ^ P^z ■> -^*) assigned 
to each EFO i. 

To select a clean DIS sample, the following cuts on global variables were applied: 

{E-Pz)tot= {E-Pz)h + E[{l-cose,) c[40,80]GeV 

y,= l-E[{l-cose,)/{2E,) <0.7 

l/jB = {E-Pz)h/{2E,) > 0.01 

Ql= {E'^sine,f/{l-y^) >2QGeV\ 

where {E-Pz)h = E.cefos ^'-Pz, 2/s = {E-Pz)h/ {E-Pz)tot [13], and 9, is the electron 
polar angle. These cuts restricted the accessible inelasticity y = Q'^/{xs) and to 
0.01 <y <0.7 and > 20 GeV^ The DIS variables x and 

were reconstructed using 

the S estimators Q| and xs = Q'^/{sy-£) [15] . 

To remove background events with isolated muons (77 /U"*"/!", J/ip and T decays) and 
residual cosmic muons, an anti-isolation cut was applied by requiring that the hadronic 
energy in a cone of radius 1 in the rj — (p plane around the muon candidate, excluding the 
muon itself, was E^^° > 0.5 GeV. From MC studies this cut was 98% (90%) efficient for 
charm (beauty). 

Jets were reconstructed from EFOs using the algorithm [H] in the longitudinally 
invariant mode About 96% of the muon candidates were associated to a jet with 
transverse momentum (including the muon) p-!^^ > 2.5 GeV and kept for further analysis. 

After the above selection, the final sample contained 11126 muons. A subsample of 35 
events with more than one muon was found, 28 of which consisted of fi~^fi~ pairs. A 
J/ip signal of 9 events was observed in the fi~^fi~ invariant mass distribution. The total 
contamination from J/i/) production was estimated with the Cascade MC, normalised 
to the observed J/ip signal. It was found to be (0.9 ± 0.3)% and was neglected in the 
analysis. 



6 Extraction of the charm and beauty fractions 

The sample of selected muon candidates contained signal muons from charm and beauty 
decays and background from in-fiight tt^ and decays and from the punch through of 



5 



hadronic jets in the muon chambers. Candidates from in-flight decays and punch through, 
which are subsequently denoted as "false muons" , were present both in the LF events and 
in events containing HQs. 

The fractions of muons originating from charm, beauty or LF events were determined 
from a simultaneous fit of three discriminating variables sensitive to different aspects of 
HQ decays: 

• the muon momentum component transverse to the axis of the associated jet, 
Pt^ = ^ P''^*I/|P'''^*|- Due to the large b mass, muons from beauty hadron decays 
have a harder p^^^ spectrum than those from charm or light quarks; 

• 6, the distance of closest approach of the muon track to the centre of the interaction 
region (beam spot) in the X, Y plane. A positive sign was assigned to S if the muon 
track crosses the axis of the associated jet in the jet hemisphere, negative otherwise. 
The beam spot position was obtained by fitting the reconstructed primary vertex 
distribution for every 2000 ep events. The size of the interaction region was 80x20 /im^ 
in Xxy. Muons from decays of long-lived heavy quarks tend to have positive 6 while 
tracks originating from the primary interaction have a symmetric 6 distribution around 
zero, corresponding to the experimental resolution. 

• p™"^*^"^, the missing transverse momentum parallel to the muon direction. The missing 
transverse momentum vector was calculated using the electron and the EFOs. The 
^missii^t ^jg^pj|-,^^jQ]2 j^g^g positlvc tall of eveuts containing semileptonic HQ decays due 
to the presence of the neutrino. 

A control sample of inclusive DIS data, selected similarly to the muon sample but without 
any muon requirement, was used to test the quality of the simulation of these variables. 
The control sample is dominated by LF events, containing, according to MC, about 
18% (1%) of c (b) events. The p^"' distribution of inclusive tracks in the control sample 
was reasonably well reproduced by both the DjangoH and the Rapgap inclusive DIS 
samples. The small differences (at most 10% at > 2 GeV) were corrected for by 
applying a bin-by-bin correction to the p^^' distribution of the LF and charm MC samples 
similarly to a previous publication [ID]. The quality of the MC description of p^^^^^^^ was 
also evaluated in the control sample by studying a similar pj^-balance variable: the missing 
transverse momentum parallel to the electron p™^''*'"''. The best description of the p™''''"° 
distribution of the inclusive DIS sample was obtained by shifting the hadronic transverse 
momentum by (0.1 ± 0.1) GeV in the MC and by increasing the hadronic transverse 
momentum resolution by (5 ± 5)% in the case of Rapgap and by (0 ± 5)% in the case 
of DjangoH. The resolution on 6 was studied using tracks in the inclusive DIS sample. 
Since it was underestimated in the MC by ^ 15%, a p-r-dependent smearing [16] was 
applied to the MC, similarly to what was done in a previous publication [17]. 
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The fractions of b, c and LF events were obtained by fitting a combination of MC dis- 
tributions to the measured three-dimensional distribution of the discriminating variables 
HH]. The fit range was \S\ < 0.1 cm, pf < 2.5 GeV and < lOGeV. A precise 

measurement of 6 was only possible inside the region covered by the BMVD. Hence for 
events with muons reconstructed in the FMUON (4% of the total) only and p™'^^"^ 
were used in the fit. A Poisson likelihood fit was used, taking into account the limited 
MC statistics. 

The global charm and beauty fractions resulting from the fit were 

fc = 0.456 ± 0.029(stat.); fb = 0.122 ± 0.013(stat.) 

with a correlation coefficient pcb = —0.43. Figure [T](a-c) shows the distributions of the 
three discriminating variables compared to the MC distributions with the normalisation 
corresponding to the fit. While S and p™'^'^"'^ provide discrimination between LF and 
HQs, discriminates between beauty and the other components. Figure 1(d) shows 
the distribution of p^"^ for a signal-enriched subsample. The distributions of p^, rj^, p'^ , 
E — Pz, Q| and for the data and for the MC samples normalised according to the fit 
are shown in Fig. [2j The overall agreement is satisfactory. 



7 Acceptance and QED corrections 

The visible cross sections for muons from charm and beauty decays, including beauty 
cascade decays via c, c, r and tp, were measured in the kinematic range 

Q2 > 20 GeV^; 0.01 <y< 0.7; pi^ > 1.5 GeV; -1.6 < 7]'' < 2.3. (1) 

The cross sections were calculated using 

where fq is the HQ fraction from the fit, is the number of reconstructed muons, Ag 
is the acceptance, Cr is the QED radiative correction, and q = c,b. Differential cross 
sections were measured by repeating the fit in bins of the reconstructed variable V as 
da/dV = a1/AVi, where is the cross section in the bin and AVi is the bin width. 

The acceptance Ag was evaluated from the MC simulation as the number of reconstructed 
muons divided by the number of true muons from decays of the quark q. This definition 
takes into account the charm and beauty events in which a "false muon" is reconstructed 
rather than a signal muon from a HQ decay. The acceptance included the efficiency of 
muon reconstruction (which in turn includes the efficiency of the muon chambers and of 
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the matching with central tracking) that was evaluated from an independent exclusive 
dimuon sample as explained in previous publications [lOl 02] . The muon reconstruction 
efficiency was around 50% for central muons with pi^ > 2 GeV. The acceptance for c (b) 
ranged from 23% (16%) at 1.5 < pf^ < 2.5 GeV to ^ 35% (25%) at pt^ > 2.5 GeV. The 
difference in acceptance between c and b was mainly due to the different contribution 
from "false muons" which was ~ 25% for c and ~ 3% for b. 

According to the MC simulation, the probability to find a "false muon" in a DIS event 
(before any muon selection) is Pfaisc ~ 0.1%, almost independently from the event being c, 
b or LF. The ability of the MC to reproduce Pfaisc was studied by comparing the number of 
LF events in the data, as obtained from the fit, to the absolute prediction by DjangoH. 
The data/MC ratio, PflfLV^faiS. was estimated as 0.80 ±0.20 in the RMUON, 1.10 ±0.20 
in the BMUON, and 1.05 ±0.40 in the FMUON, in agreement with previous studies |49j . 

The cross sections were corrected to the QED Born level, calculated using a running 
coupling costant aem, such that they can be compared directly to the QCD predictions 
from the HvQDis program (Section [2]). The radiative corrections were obtained as Cr = 
o"Born/o"rad7 wherc o"Born IS the Rapgap cross sectiou with the QED corrections turned 
off but keeping acm running and cxrad is the Rapgap cross section with the full QED 
corrections, as in the standard MC samples. The corrections were typically Cr ~ 1.05 
and at maximum 1.10 in the highest bin. 



8 Systematic uncertainties 

The following systematic uncertainties were considered (the effects on the total visible 
cross section for c and for b is given in parentheses): 

1. B/RMUON efficiency: it was varied by its uncertainty of on average ±5% (=f5, ±5)%; 

2. FMUON efficiency: it was varied by ±20% (^2, ^5)%; 

3. "false muon" probability: it was varied within the corresponding uncertainty for each 
muon detector (],]4,=f1)%; 

4. global energy scale: it was varied by ±2% (^5,^2 )%j 

5. calibration of p™^''*'"^: it was evaluated by varying the hadronic transverse momentum 
in the MC by ±0.1 GeV, as allowed by the transverse momentum balance in the 
control sample (±12,^1)%; 

6. hadronic energy resolution: it was varied in the MC by ±5% as allowed by the trans- 
verse momentum balance in the control sample ±7)%; 
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7. simulation of the tails of p™'^^"'^: the fits were redone in the restricted range < 
5GeV (0,-6)%; 

8. resolution on 6: the smearing applied to the MC was varied by ±25% as allowed by 
the control sample (+25-9^ 

9. shape of LF and charm: it was evaluated by varying the p^^' correction by ±50% 
(±1.5+^)%; 

10. hadronic energy flow near the muon: it was evaluated by varying the cut on E^^° by 
t°;iGeV(0,-o^)%; 

11. jet description: the cut on p!^^ was varied by ±0.5 GeV (±2.5,^2^5 

12. charm SL decay spectrum: the reweighting to the CLEO model was varied by ±50%, 

(;3,l2)%; 

13. MC model dependence: the Rapgap c and b samples were reweighted to reproduce 
the corresponding measured differential cross sections in or in pt^ and the largest 
deviation from the nominal cross section was taken (±6, ±20)%; 

14. higher order effects: this uncertainty was evaluated by varying the HQ distribution 
before parton showering in Rapgap by the difference between NLO and leading order, 
as evaluated with HvQDis {tto,tl )%; 

15. MVD efficiency: the efficiency of the cut on the number of MVD hits of (90 ± 3)% 
was varied by its uncertainty (=f3, ±3)%; 

16. CTD simulation: tracks were required to pass > 4 superlayers in the B/RMUON 
region and to have > 7 degress of freedom in the FMUON region (±1, 0)%; 

17. integrated luminosity: measurement uncertainty (=f2. 6, =f2. 6)%. 

The above uncertainties were summed in quadrature to obtain the total systematic un- 
certainty 

9 Cross sections 

The visible cross sections for muons from charm and beauty decays in the kinematic region 
of Eq. ([T]) are 

= 164 ± 10(stat.) +3?(syst.) pb 
a'' = 63 ± 7(stat.) +i^(syst.) pb 

to be compared with the NLO QCD cross sections obtained with HvQDiS of a'^ = 
184^*^40 pb and = 33 ± 5 pb. The agreement is good for charm while the beauty cross 
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section is 2.3 (1.9) standard deviations above the central (upper) HvQDis result. The 
visible cross sections are a factor 1.04 and 2.27 higher than the Rapgap MC predictions 
for c and b, respectively. 

The differential cross sections as a function of p^, r/^, Q^, and x are presented in Tabled] 
and compared in Fig. [3] to the NLO QCD predictions based on HvQDis. The Rapgap 
MC predictions are also shown, normalised according to the result of the global fit. Charm 
and beauty cross sections are similar for pt^ > 3.5 GeV. 

The charm cross sections are in good agreement with the HvQDis calculations. The ten- 
dency of the beauty cross section to lie above the central NLO prediction is concentrated 
at low pi^ and low Q^. The statistical significance of the difference between the data and 
the NLO predicitons is similar to that obtained for the total visible cross section since the 
uncertainties are dominated by correlated systematics. 

Both NLO calculations and the Rapgap MC give in general a good description of the 
shape of the differential cross sections. The distributions for beauty is somewhat 
steeper than predicited by Rapgap and HvQDis. 



10 Extraction of 

The heavy quark contribution to the proton structure functions, Fg^'^, and the reduced 
cross section a'^^ are defined in analogy with the inclusive case from the double differential 
cross section in x and for the production of the quark q: 



Ff{x,Q')-^Ff{x,Q') 



/C a'"^{x,Q\s 



dx dQ"^ 

where K. = Y+{27ra^^J/{xQ^) and r+ = 1 + (1 - y)\ 

The muon cross sections, o"^, measured in bins of x and Q^, were used to extract at 
a reference points in the x, plane by: 

^''''''ix,Q^) 



Ff(x,Q2) = a^ 



where F2^^'*''(x, Q^) and a^'*'^ were calculated at NLO in the FFNS using the HvQDis 
program. The reference points were chosen close to the average x and of the events 
within each bin. Charm produced in b decays was not included in the definition of cr'^. 

This procedure contains several corrections: the extrapolation from the restricted muon 
kinematic range {pf^ > l.SGeV, — 1.6 < t]'^ < 2.3) to the full muon phase space; the 
g — >■ /i branching ratio; the correction for the longitudinal structure function Fl'^ and the 
correction from a bin-averaged cross section to a point value (bin centring). 
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The largest uncertainty is related to the extrapolation to the full muon phase space. The 
kinematic acceptance, A, defined as the fraction of muons from HQ decays that was 
generated in the restricted kinematic region is, on average, (.4) = 13%(27%), for charm 
(beauty). According to HvQDis, in the charm case, A becomes sizeable {A > 0.25(A)) 
when one of the two charm quarks in the event has p-p > 3 GeV and its rapidity is in 
the range (—1.5 : 2.5), which corresponds to the phase space containing 88% of the cross 
section. In the beauty case, A is sizeable over the full HQ phase space. 

The theoretical uncertainty in the extraction of ^2^'^ was evaluated by varying the HvQDis 
parameters as explained in Section [2] and by using a different PDF set (CTEQ5F). It is 
dominated by the fragmentation uncertainty. As a further check, was also evaluated 
taking A from Rapgap and found to be consistent within the quoted uncertainties. 

The muon cross sections in bins of x and are given in Table [21 The extracted F2'^ 
and are presented in Tables [3] and H] and shown in Figs. H] and O Also given in 
Tables [3] and H] are the factor A and the correction for the longitudinal structure function 
Cl = a'^^ / F^'^ as obtained from the NLO theory. The effects of the individual sources of 
systematic and theoretical uncertainty are given in Tables [5] and [6] in the Appendix [^. 

Figure m also contains a comparison of Fg^^ with previous results based on the measurement 
of D* mesons from ZEUS [2] and to results from the HI collaboration based on inclusive 
hfetime tagging (VTX) [3]. The previous results were corrected to the values of the 
present analysis, using the NLO theory. The agreement of the different data sets, obtained 
with different charm tagging techniques, is good. At high Q^, the precision of present 
data is similar or better than for the previous results. The NLO QCD calculations are 
also shown. 

Figure E] shows the extracted F!^^ from this analysis and also a previous HI result [3], 
corrected to the reference values used in the present analysis. The two data sets are 
in good agreement. The precision of the present measurement is similar to that of the HI 
data at high Q^. The QCD calculations are also shown. 

The structure functions -F^^^ and F2^ are also presented in Figs. [6] and [7] as functions of 
for fixed values of x, compared to previous results corrected to the same reference x 
used in the present analysis. 



^ The effects of the individual sources of systematic and theoretical uncertainty are also available from 
http : //www-zeus . desy . de/publicjresults/f unctiondb.php?id=ZEUS-pub-09-003. 
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11 Summary 



The production of charm and beauty quarks was measured in DIS using their decay into 
muons. Total and differential cross sections for muons from c and b decays were measured 
in the kinematic region 

Q2 > 20 GeV^; 0.01 <y< 0.7; > 1.5 GeV; -1.6 < < 2.3 

and compared to NLO QCD calculations. The agreement is good for charm. Beauty is 
about a factor two above the central QCD prediction although still compatible within 
statistical and systematic uncertainties. The heavy quark contribution to the proton 
structure function F2 was also measured and found to agree well with other measurements 
based on independent techniques. For > 60 GeV^ the present results are of comparable 
or higher precision than those previously existing. 



12 Acknowledgements 

We appreciate the contributions to the construction and maintenance of the ZEUS de- 
tector of many people who are not listed as authors. The HERA machine group and the 
DESY computing staff are especially acknowledged for their success in providing excel- 
lent operation of the collider and the data analysis environment. We thank the DESY 
directorate for their strong support and encouragement. It is also a pleasure to thank 
Pedro Jimenez Delgado, Stefano Forte, Albero Guffanti, Eric Laenen, Pavel Nadolsky, 
Jack Smith and Paul Thompson for many illuminating discussions. 



12 



References 



[1] ZEUS Coll., J. Breitweg et al, Phys. Lett. B 407, 402 (1997); 
ZEUS Coll., J. Breitweg et al, Eur. Phys. J. C 12, 35 (2000); 
HI Coll., C. Adloff et al, Phys. Lett. B 528, 199 (2002); 
HI Coll., A. Aktas et al, Eur. Phys. J. C 38, 447 (2005). 

ZEUS Coll., S. Chekanov et al., Phys. Rev. D 69, 012004 (2004). 

HI Coll., A. Aktas et al, Eur. Phys. J. C 40, 349 (2005); 
HI Coll., A. Aktas et al, Eur. Phys. J. C 45, 23 (2006). 

ZEUS Coll., S. Chekanov et al., Phys. Lett. B 599, 173 (2004). 

HI Coll., A. Aktas et al., Eur. Phys. J. C 41, 453 (2005). 

J.D. Bjorken, Phys. Rev. 179, 1547 (1969). 

J. Smith and W.L. van Neerven, Nucl. Phys. B 374, 36 (1992). 

B.W. Harris and J. Smith, Nucl. Phys. B 452, 109 (1995). 

B. W. Harris and J. Smith, Phys. Lett. B 353, 535 (1995). Erratum-ibid B 359 
(1995) 423. 

ZEUS Coll., S. Chekanov et al., Phys. Rev. D 67, 012007 (2003). 

C. Peterson et al., Phys. Rev. D 27, 105 (1983). 
P. Nason and C. Oleari, Nucl. Phys. B 565, 245 (2000). 

M. Cacciari, P. Nason and C. Oleari, JHEP 0604, 006 (2006). Oct 2005. 48p. 

BELLE Coll., R. Seuster et al., Phys. Rev. D 73, 032002 (2006). 

CLEO Coll., N. E. Adam et al., Phys. Rev. Lett. 97, 9 (2006). 

T. Sjostrand et al., Comp. Phys. Comm. 135, 238 (2001). 

Particle Data Group, W.-M. Yao et al., J. Phys. G 33, 1 (2006). 

BELLE Coll., K. Abe et al., Phys. Lett. B 547, 181 (2002); 
BABAR Coll., B. Aubert et al., Phys. Rev. D 67, 031101 (2002). 

ZEUS Coll., S. Chekanov et al., JHEP 0904, 082 (2009); 
HI Coll., F.D. Aaron et al., Eur. Phys. J. C 59, 589 (2009). 

M. Botje, (^Ci^A'^f/M werszon (unpublished), available on 

http : //www.nikhef .nl/~h24/qcdcode/index . html. 

H. Jung, Comp. Phys. Comm. 86, 147 (1995). 



13 



[22] A. Kwiatkowski, H. Spiesberger and H.-J. Mohring, Comp. Phys. Comm. 

69, 155 (1992). Also in Proc. Workshop Physics at HERA, eds. W. Buchmiiller and 

G. Ingelman, (DESY, Hamburg, 1991). 

[23] CTEQ Coll., H.L. Lai et al., Eur. Phys. J. C 12, 375 (2000). 

[24] G.A. Schuler and H. Spiesberger, Proc. Workshop on Physics at HERA, 

W. Buchmiiller and G. Ingelman (eds.). Vol. 3, p. 1419. Hamburg, Germany, DESY 
(1991); 

H. Spiesberger, HERACLES and djangoh.- Event Generation for ep Interactions at 
HERA Including Radiative Processes, 1998, available on 

http : / / www . desy . de/~hspiesb/ dj angoh . html. 

[25] G. Ingelman, A. Edin and J. Rathsman, Comp. Phys. Comm. 101, 108 (1997). 

[26] H. Jung and G.P. Salam, Eur. Phys. J. C 19 (2001). 

[27] ZEUS Coll., S. Chekanov et al., Eur. Phys. J. C 44, 13 (2005). 

[28] T. Sjostrand, Comp. Phys. Comm. 82, 74 (1994). 

[29] R. Brun et al., geantS, Technical Report CERN-DD/EE/84-1, CERN, 1987. 

[30] ZEUS Coll., U. Holm (ed.). The ZEUS Detector. Status Report (unpublished), 

DESY (1993), available on http://www-zeus.desy.de/bluebook/bluebook.html. 

[31] A. Pohni et al, Nucl. Inst. Meth. A 581, 31 (2007). 

[32] N. Harnew et al, Nucl. Inst. Meth. A 279, 290 (1989); 

B. Foster et al., Nucl. Phys. Proc. Suppl. B 32, 181 (1993); 
B. Foster et al., Nucl. Inst. Meth. A 338, 254 (1994). 

[33] M. Derrick et al., Nucl. Inst. Meth. A 309, 77 (1991); 
A. Andresen et al, Nucl. Inst. Meth. A 309, 101 (1991); 
A. Caldwell et al., Nucl. Inst. Meth. A 321, 356 (1992); 
A. Bernstein et al., Nucl. Inst. Meth. A 336, 23 (1993). 

[34] A. Bamberger et al., Nucl. Inst. Meth. A 382, 419 (1996). 

[35] A. Bamberger et al., Nucl. Inst. Meth. A 401, 63 (1997). 

[36] G. Abbiendi et al, Nucl. Inst. Meth. A 333, 342 (1993). 

[37] J. Andruszkow et al.. Preprint DESY-92-066, DESY, 1992; 
ZEUS Coll., M. Derrick et al, Z. Phys. C 63, 391 (1994); 
J. Andruszkow et al. Acta Phys. Pol. B 32, 2025 (2001). 

[38] M. Helbich et al, Nucl. Inst. Meth. A 565, 572 (2006). 

[39] P.D. Allfrey et al., Nucl. Inst. Meth. A 580, 1257 (2007); 

W. H. Smith, K. Tokushuku and L. W. Wiggers, Proc. Computing in High-Energy 



14 



Physics (CHEP), Annecy, France, Sept. 1992, C. Verkerk and W. Wojcik (eds.), 
p. 222. CERN, Geneva, Switzerland (1992). Also in preprint DESY 92-150B. 

[40] ZEUS Coll., S. Chekanov et al., Phys. Rev. D 70, 12008 (2004). 

[41] G.M. Briskin, Diffractive Dissociation in ep Deep Inelastic Scattering. Ph.D. 
Thesis, Tel Aviv University, Report DESY-THESIS 1998-036, 1998. 

[42] M. Turcato, Measurement of Beauty Photoproduction at HERA. Ph.D. Thesis, 
Universita di Padova, Report DESY-THESIS-03-039, 2002. 

[43] U. Bassler and G. Bernardi, Nucl. Inst. Meth. A 361, 197 (1995). 

[44] S. Catani and M.H. Seymour, Nucl. Phys. B 485, 291 (1997). Erratum in 
Nucl. Phys B 510, 503 (1997). 

[45] S. Catani et al., Nucl. Phys. B 406, 187 (1993); 

S.D. Ellis and D.E. Soper, Phys. Rev. D 48, 3160 (1993). 

[46] S. Miglioranzi, Beauty Photoproduction at HERA U with the ZEUS Experiment. 
PhD Thesis, University College London, 2006, available on 

http : //www-zeus . desy.de/physics/hfla/public/Theses/PhD_niiglioranzi .pdf . 
[47] ZEUS Coll., S. Chekanov et al., JHEP 0904, 133 (2009). 

[48] M. Bindi, Measurement of the Charm Production Cross Section in DIS Events at 
HERA. PhD Thesis, Universita di Bologna, Report DESY-THESIS-2008-026, 2008. 

[49] A. Longhin, Measurement of Beauty Production at HERA with a D* + fi Tag. 
Ph.D. Thesis, Universita di Padova, Report DESY-THESIS-2004-050, 2003. 



15 



Pt 


da'^/ dpt^ 


A Stat. 


^syst. 


da^/ dpt^ 


^stat. 


^syst. 


Pc,b 


(GeV) 




(pb/GeV) 






(pb/GeV) 






1.5 : 2.5 


113 


±10 


+21 

-23 


48 


±8 


+15 
-13 


-0.48 


2.5 : 3.5 


32.8 


±3.8 


+5.7 
-5.4 


15.4 


±2.7 


+4.1 
-2.8 


-0.46 


3.5 : 5.0 


6.0 


±1.4 


+1.4 
— 1.6 


6.3 


±1.1 


+0.9 
—0.7 


-0.48 


5.0 : 10.0 


0.97 


±0.21 


+0.13 
-0.20 


0.76 


±0.22 


+0.10 
-0.15 


-0.45 




da^dr]'' 


^stat. 
(pb) 


^syst. 


da''/d7]f' 


^stat. 
(pb) 


^syst. 


Pc,b 


-1.60 : -0.90 


20.4 


±3.3 


+4.2 
—4.1 


5.0 


±2.0 


+1.9 
— 1.6 


-0.36 


-0.90 : -0.40 


40.7 


±6.0 


+7.7 
—8.0 


13.6 


±3.5 


+4.1 
— 1.9 


-0.20 


-0.40 : +0.00 


60.9 


±7.8 


+11.8 
— 14.2 


17.3 


±4.8 


+4.9 
—3.5 


-0.37 


+0.00 : +0.50 


67.0 


±7.1 


+10.2 
— 12.3 


21.2 


±4.4 


+6.4 
—3.6 


-0.42 


+0.50 : +1.48 


47.7 


±6.4 


+8.9 
-8.8 


20.7 


±4.3 


+6.0 
-5.5 


-0.49 


+ 1.48 : +2.30 


33.4 


±10.0 


+15.3 
-8.9 


16.4 


±6.9 


+5.7 
-8.4 


-0.41 




da^dQ^ 


^stat. 


^syst. 


da''/dQ^ 


^stat. 


^syst. 


Pc,b 


(GeV^) 




(pb/GeV^) 






(pb/GeV2) 






20 : 40 


3.43 


±0.40 


+0.72 
—0.66 


1.46 


±0.24 


+0.30 
—0.39 


-0.44 


40 : 80 


1.22 


±0.13 


+0.18 
—0.22 


0.546 


±0.086 


+0.109 
—0.098 


-0.41 


80 : 200 


0.289 


±0.031 


+0.053 
-0.054 


0.124 


±0.023 


+0.020 
-0.019 


-0.36 


200 : 500 


0.0447 


±0.0071 


+0.0050 
-0.0083 


0.0131 


±0.0049 


+0.0035 
-0.0024 


-0.47 


500 : 10000 


0.00063 


±0.00013 


+0.00014 
-0.00010 


0.00013 


±0.00008 


+0.00005 
-0.00003 


-0.38 


X 


c/cr^/ dx 


^stat. 

(nb) 


^syst. 


da''/ dx 


^stat. 

(nb) 


^syst. 


Pc,b 


0.0003 : 0.0010 


35.3 


±5.6 


+10.2 
-6.2 


17.4 


±3.9 


+3.5 
-3.6 


-0.16 


0.0010 : 0.0020 


35.2 


±4.1 


+4.6 
-7.1 


12.4 


±2.6 


+2.7 

-3.0 


-0.39 


0.0020 : 0.0040 


16.1 


±2.2 


+3.8 
-3.7 


8.0 


±1.4 


+1.6 
-1.2 


-0.51 


0.0040 : 0.0100 


7.38 


±0.72 


+1.22 
-1.30 


2.04 


±0.45 


+0.52 
-0.39 


-0.45 


0.0100 : 0.1000 


0.417 


±0.050 


+0.068 
-0.081 


0.076 


±0.028 


+0.030 
-0.014 


-0.53 



Table 1: Muon differential cross sections for charm and beauty as a function 
of T]^, pi^p, and x. The last column shows the statistical correlation coefficient 
between charm and beauty. 



16 



bin 






X 






Astat. 


A syst. 




Astat. 


A syst. 


Pc,b 




(GeV2) 








(pb) 






(pb) 






1 


20 


: 60 


0.0003 : 


0.0012 


32.9 


±4.6 


+ 8.1 

-5.9 


13.9 


±2.9 


+3.2 
-2.3 


-0.29 


2 


20 


: 60 


0.0012 : 


0.0020 


17.7 


±3.1 


+1.9 
-4.9 


5.7 


±2.0 


+1.4 
-1.2 


-0.42 


3 


20 


: 60 


0.0020 : 


0.0035 


16.2 


±3.3 


+3.7 
-3.0 


5.5 


±2.0 


+1.8 
-1.6 


-0.51 


4 


20 


: 60 


0.0035 : 


0.0060 


35.1 


±5.7 


+10.9 
-7.2 


7.9 


±3.6 


+4.2 
-4.0 


-0.56 


5 


60 


: 400 


0.0009 : 


0.0035 


17.2 


±2.7 


+3.8 
-2.7 


8.8 


±1.9 


+1.6 
-1.4 


-0.38 


6 


60 


: 400 


0.0035 : 


0.0070 


18.4 


±2.3 


+3.0 
-3.4 


4.2 


±1.5 


+1.2 
-0.9 


-0.35 


7 


60 


: 400 


0.0070 : 


0.0400 


33.6 


±3.5 


+6.1 
-6.4 


8.6 


±2.3 


+2.1 
-2.1 


-0.46 


8 


400 


:10000 


0.0050 : 


1.0000 


7.6 


±1.5 


+1.2 
-1.4 


1.6 


±0.9 


+0.4 
-0.4 


-0.45 



Table 2: Muon cross sections for charm and beauty in bins of and x. The 
last column shows the statistical correlation coefficient between charm and beauty. 



bin 


(GeV2) 


X 


T7CC 

^2 


^stat. 


^syst. 


^theo. 


A 


Cl 


1 


30 


0.0008 


0.318 


±0.044 


+0.078 
-0.057 


+0.061 
-0.042 


0.096 


0.980 


2 


30 


0.0016 


0.219 


±0.038 


+0.024 
-0.061 


+0.043 
-0.016 


0.114 


0.996 


3 


30 


0.0025 


0.176 


±0.036 


+0.040 
-0.033 


+0.032 
-0.021 


0.113 


0.998 


4 


30 


0.0055 


0.143 


±0.023 


+0.044 
-0.029 


+0.028 
-0.009 


0.096 


1.000 


5 


130 


0.0025 


0.298 


±0.047 


+0.066 
-0.046 


+0.044 
-0.025 


0.175 


0.955 


6 


130 


0.0055 


0.228 


±0.029 


+0.037 
-0.042 


+0.030 
-0.015 


0.220 


0.993 


7 


130 


0.0130 


0.151 


±0.016 


+0.027 
-0.029 


+0.021 
-0.011 


0.209 


0.999 


8 


1000 


0.0300 


0.114 


±0.023 


+0.018 
-0.021 


+0.010 
-0.007 


0.371 


0.987 



Table 3: The structure function F2'^{x, Q^). The last two columns show the muon 
kinematic acceptance, A, and the correction for the longitudinal structure function, 
Cl. 
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bin 




X 




^stat. 


^syst. 


^thco. 


A 


Cl 


1 


30 


0.0008 


0.0220 


±0.0047 


+0.0049 
-0.0037 


+0.0011 
-0.0010 


0.260 


0.992 


2 


30 


0.0016 


0.0131 


±0.0047 


+0.0032 
-0.0028 


+0.0009 
-0.0003 


0.264 


0.998 


3 


30 


0.0025 


0.0114 


±0.0043 


+0.0037 
-0.0034 


+0.0005 
-0.0004 


0.251 


0.999 


4 


30 


0.0055 


0.0080 


±0.0036 


+0.0041 
-0.0041 


+0.0004 
-0.0003 


0.189 


1.000 


5 


130 


0.0025 


0.0489 


±0.0105 


+0.0088 
-0.0076 


+0.0024 
-0.0018 


0.300 


0.962 


6 


130 


0.0055 


0.0175 


±0.0064 


+0.0052 
-0.0039 


+0.0007 
-0.0007 


0.319 


0.994 


7 


130 


0.0130 


0.0149 


±0.0039 


+0.0037 
-0.0037 


+0.0007 
-0.0006 


0.281 


0.999 


8 


1000 


0.0300 


0.0104 


±0.0061 


+0.0028 
-0.0025 


+0.0004 
-0.0004 


0.420 


0.983 



Table 4: The structure function F2^{x, Q^). The last two columns show the muon 
kinematic acceptance, A, and the correction for the longitudinal structure function, 
Cl- 



18 



ZEUS 




■10 -5 5 10 

pmissim ^Qg^j 



(0 

c 
o 

^10^ 


- 1 1 1 1 1 1 1 1 

^ (b) r 


1 1 1 1 1 1 1 1 1 1 1 - 

ri - 

f "i 


10^ 


• 1 -! 

If ' 


[■■ •^^^^^^^^^^^^ 


10 


r#J' ■- 1"' --. 

=-F-l""i 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 i- 



-0.1 -0.05 0.05 0.1 

6 (cm) 




(0 

c 
o 



T 1 1 1 1 1 1 1 1 r 

FMUON or 5>0.01 cm '. 
pmissim^2 GeV 




pT' (GeV) 



Pt (GeV) 



Figure 1: Distributions of (a) , (b) 6, (c) p^'f^^ for the selected sample of 

muons in DIS, and of (d) p^?' for a signal- enriched subsample with p^^^^^^^ > 2 GeV 
and either a muon in FMUON or 6 > 0.01 cm. The data (points) are compared 
to the MC expectation (solid line) with the normalisation of the c (dotted line), b 
(shaded histogram) and light flavours (dashed line), LF, components obtained from 
the global fit. The error bars correspond to the square root of the number of entries. 
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Figure 2: Distributions of (a) p!^, (h) rf^ , (c) pJ^*, (d) {E — Pz)tot, (^) Q| 
(f) ^T. for the selected sample of muons in DIS. The data (points) are compared 
to the MC expectation with the normalisation of the c, h and light flavours, LF, 
components obtained from the global fit. 
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Figure 3: Differential muon cross sections for c and h as a function of (a) pf^, 
(b) ?7^, (c) Q^, and (d) x. The inner error bars show the statistical uncertainty 
while the outer error bars show the systematic and statistical uncertainties added in 
quadrature. The bands show the NLO QCD predictions obtained with the HvQDis 
program and the corresponding uncertainties. The differential cross sections from 
Rapgap, scaled by the factors corresponding to the result of the global fit (1.04 for 
c and 2.27 for h), are also shown. 
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Figure 4: The structure function F^'^ (filled symbols) compared to previous results 
(open symbols) and to the NLO QCD predictions in the FFNS using the ZEUS-S 
PDF fit. The inner error bars are the statistical uncertainty while the outer bars 
represent the statistical, systematic and theoretical uncertainties added in quadra- 
ture. The band represents the uncertainty on the NLO QCD prediction. Previous 
data have been corrected to the reference values used in this analysis: ZEUS D* 
500 -> 1000 GeV^; HI VTX 25 ^ 30, 200 130, 650 ^ 1000 GeV'^ . 
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Figure 5: The structure function F!^^ (filled symbols) compared to previous results 
(open symbols) and to the NLO QCD predictions in the FFNS using the ZEUS-S 
PDF fit. The inner error bars are the statistical uncertainty while the outer bars 
represent the statistical, systematic and theoretical uncertainties added in quadra- 
ture. The band represents the uncertainty on the NLO QCD prediction. Previous 
data have been corrected to the reference values used in this analysis: 25 30, 
200 ^ 130, 650 1000 GeV^. 
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ZEUS 




Figure 6: The structure function Fg'^ (filled symbols ) plotted as a function ofQ^ for 
fixed X values. The curves represent the NLO QCD predictions in the FFNS using 
the ZEUS-S PDF fit. The inner error bars are the statistical uncertainty while the 
outer bars represent the statistical, systematic and theoretical uncertainties added 
in quadrature. The band represents the uncertainty on the NLO QCD prediction. A 
selection of previous data (open symbols) is also shown, corrected to the reference 
X values used m this analysis: ZEUS D* : 0.001 0.0008, 0.0015 0.0016, 
0.003 0.0025, 0.006 0.0055, 0.012 0.013,- HI VTX: 0.0005 0.0008, 
0.002 ^ 0.0025, 0.005 ^ 0.0055, 0.032 ^ 0.030. 
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ZEUS 




Figure 7: The structure function F!^'^ (filled symbols) plotted as a function ofQ^ for 
fixed X values. The curves represent the NLO QCD predictions in the FFNS using 
the ZEUS-S PDF fit. The inner error bars are the statistical uncertainty while the 
outer bars represent the statistical, systematic and theoretical uncertainties added 
in quadrature. The band represents the uncertainty on the NLO QCD prediction. 
All the previous data (open symbols) are also shown, corrected to the reference x 
values used m this analysis: 0.0005 0.0008, 0.002 ^ 0.0025, 0.005 ^ 0.0055, 
0.032 ^ 0.030. 
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Appendix: 

Tables of systematic and theoretical uncertainties 



Syst. 








Ff 


bin 














pbb 


bin 










1 


2 


3 


4 


5 


6 


7 


8 


1 


2 


3 


4 


5 


6 


7 


8 


la 


-5 


-6 


-6 


-4 


-6 


-6 


-6 


-6 


-6 


-5 


-2 


7 


-6 


-6 


-2 


-5 


lb 


6 


7 


6 


4 


6 


6 


6 


7 


7 


5 


1 


-7 


6 


6 


1 


4 


2a 





1 


-1 


-6 


1 





-1 


1 


-2 


-6 


-12 


-37 


-2 


-4 


-14 


-6 


2b 





-1 


1 


5 


-1 





1 


-1 


2 


5 


11 


43 


2 


3 


14 


6 


3a 


-3 


-2 


-3 


-3 


-4 


-3 


-2 


-2 





-2 





-8 


1 


1 


-2 


-3 


3b 


3 


2 


3 


4 


4 


3 


2 


1 





3 


1 


11 


-1 


-1 


4 


5 


4a 


2 


-11 


8 


-3 


-1 


-6 


-6 


-9 


2 


7 


-8 


14 


-7 


10 


1 


-8 


4b 


-2 


-1 


2 


5 


15 


2 


5 


12 


9 


-9 


-8 


13 


-7 


11 





-4 


5a 


4 


4 


17 


20 


9 


11 


12 


3 


2 


-3 





5 


-2 


4 


-6 


3 


5b 


-4 


-20 


-10 


-15 


-1 


-9 


-11 


-2 


5 


11 


-1 


10 


-3 


-9 


3 


-5 


6a 


19 


-5 


1 


1 





-2 


-10 





3 


5 


-13 


-28 


-2 


-9 


-3 


-7 


6b 


-9 





3 


19 


4 


2 


8 


3 


6 


-3 


22 


3 


2 


11 


-1 


4 


7 


8 


-5 


-2 


1 


-5 


-4 


4 


4 


1 


-8 


7 


-8 


-4 


11 


-16 


-9 


8a 





-2 


-1 


1 





2 


-1 


-3 


11 


6 


-4 


-1 


8 


6 


11 


21 


8b 





4 


-3 


3 


3 





-3 


4 


-11 


-10 


-5 


-10 


-5 


-7 


-1 


-16 


9a 


2 


-1 


-4 


1 





-1 


-2 


-1 


-5 


7 


14 


2 


5 


10 


8 


6 


9b 


-1 


1 


3 


-1 





1 


1 





4 


-4 


-9 


-1 


-3 


-6 


-5 


-3 


10a 


-2 


4 


1 


1 


-2 


-1 








-1 


-4 


-1 


-1 


-1 


-1 


2 


1 


10b 


1 


-5 





2 


-2 


-1 


2 


-1 


-4 


3 


4 


-5 


6 


-6 





5 


11a 





-1 


3 


9 


6 


-2 


1 





-5 


6 


11 


-8 


-3 


-2 


-1 


3 


lib 





-7 


-9 


-4 


2 


2 








1 


5 


3 


6 


-2 


-3 


1 


2 


12a 


-3 





-2 





-3 


-3 


-2 





2 


-2 


2 


-7 


4 


6 


4 


5 


12b 


2 





1 


-1 


2 


2 


2 


-1 


-1 


2 


-1 


9 


-3 


-5 


-3 


-3 


13 


-1 


-1 


-2 


1 


-2 





-1 





9 


10 


5 


-2 


9 


5 


6 


1 


14a 


6 


1 


3 


1 


4 


6 


3 





4 


5 


5 


3 


-2 


-8 


6 


8 


14b 


-12 


-7 


-8 


-9 


-11 


-12 


-7 


-10 


-6 


-8 


-11 


4 


3 


12 


-6 


2 


15a 


-3 


-3 


-3 


-2 


-3 


-3 


-3 


-3 


-3 


-2 


-1 


3 


-3 


-2 


-1 


-2 


15b 


3 


3 


3 


2 


3 


3 


3 


3 


3 


2 


1 


-3 


3 


2 


1 


2 


16 


6 


-2 


5 


-5 


4 


3 


-2 


-8 


1 


6 


-12 


-5 





-6 


5 


2 


17a 


-3 


-3 


-3 


-3 


-3 


-3 


-3 


-3 


-3 


-3 


-3 


-3 


-3 


-3 


-3 


-3 


17b 


3 


3 


3 


3 


3 


3 


3 


3 


3 


3 


3 


3 


3 


3 


3 


3 



Table 5: Systematic uncertainties of the and measurements. The first 
column gives the systematic variation number as reported in Section 8, with "a" 
and "b" corresponding to variations in opposite directions. The other columns list 
the effect of each variation on the measurements in percent. 
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Syst. 










bin 
















bin 










1 


2 


3 


4 


5 


6 


7 


8 


1 


2 


3 


4 


5 


6 


7 


8 




4 


3 


1 


2 


1 


1 








1 


2 





-1 








-1 





/2 


-8 


-2 


-7 





-4 


-2 





1 


-1 





-1 





-1 





1 







1 


-1 


-3 


-2 


-2 


-2 


-3 


-2 





1 





-2 








-1 


-1 


/2 


-1 





2 


6 


1 


2 


5 


3 


-1 


1 


1 


1 


1 





1 


1 


mg+ 


-8 


-5 


-3 


-2 


-3 


-1 


-1 


-2 


-2 


-1 


-1 





-1 


-1 





-1 




5 


4 


3 


3 


3 


2 


2 


1 


2 


2 


1 





1 


1 


2 





PDF + 


2 


-1 


-1 





1 





2 








1 


-1 


-1 








-1 


-1 




-1 


1 


-2 


-1 


-1 














1 





1 








1 





CTEQ 


4 


1 


-1 


1 


1 





2 


2 





1 


-1 


3 





1 


2 


1 


e+ 


16 


17 


17 


17 


13 


12 


12 


7 


2 


4 


2 


2 


3 


3 


3 


2 




-6 


-2 


-8 


-4 


-5 


-4 


-5 


-4 


-3 


-1 


-3 


-2 


-3 


-3 


-3 


-2 


S + PII 


5 


6 


4 


4 


3 


3 


2 





3 


5 


3 


2 


3 


3 


2 







-4 


-4 


-4 


-4 


-4 


-4 


-4 


-4 


-2 


-2 


-2 


-2 


-2 


-2 


-2 


-2 




4 


4 


4 


4 


4 


4 


4 


4 


2 


2 


2 


2 


2 


2 


2 


2 



Table 6: Theoretical uncertainties of the Fg'^ and measurements. The first 
column gives the parameter varied in the calculation as reported in Sections 2 and 
10: the factorisation (fip) o,nd renormalisation (fiR,) scales, the HQ mass (mq), 
the variation of the ZEUS PDF by its uncertainty (PDF), the use of the CTEQ5 
PDF (CTEQ), the Peterson fragmentation parameter (e), the use an alternative 
fragmentation variable (E+p\\) and the SL branching ratio (B). The other columns 
list the effect of each variation on the measurements in percent. 
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